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In this work, we prove a generalization of Osgood's test for the ex- 
plosion of the solutions of initial-value problems. We also establish a 
comparison criterion for the solution of integral equations with noise, and 
provide estimations of the time of explosion of problems arising in the 
investigation of crack failures where the noise is the absolute value of the 
Brownian motion. 
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1 Introduction 

Let xq be a positive, real number, let 6 be a positive, real- valued function defined 
on [0, oo), and suppose that y is an extended- real- valued function with the same 
domain as b. The present work is motivated by a criterion for the explosion of 
the solutions of ordinary differential equations of the form 
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Abstract 




t > 0, 



(1) 



More precisely, the time of explosion of the solution of this initial- value problem 
is the nonnegative, extended-real number = sup{t > : y(t) < oo}. The 
above-mentioned criterion is called Osgood's test after its author [Kl, and it 
states that t,, is finite if and only if J ds/b{s) < oo. In such case, = 

Cods/Hs). 

A natural question readily arises about the possibility to extend Osgood's 
test to more general, initial- value problems, say, to problems in which the drift 
function b in the ordinary differential equation of ([T]) is multiplied by a suit- 
able, nonnegative function of t. Another direction of investigation would be to 
investigate conditions under which the solutions of the integral form of such 
equation with a noise function added, explode in a finite time (see for example 
[9]). Evidently the consideration of these two problems as a single one is an 
interesting topic of study per se. In fact, the purpose of this paper is to provide 
a generalization of Osgood's test to integral equations with noise, which gener- 
alize the problem presented in ([T]). Important, as it is in the recent literature 
[3 ini [Zl E] , the problem of establishing analytical conditions under which the 
time of explosion of the problem under investigation, is tackled here. In the 
way, we establish a comparison criterion for the solutions of integral equations 
with noise, and show some applications to the spread of cracks in rigid surfaces. 

Our manuscript is divided in the following way: Section [2] introduces the 
integral equation with noise that motivates this manuscript, along with a con- 
venient simplification for its study; a generalization of Osgood's test is presented 
in this stage for the associated initial- value problem for both scenarios: noiseless 
and noisy systems. Section [3] establishes a comparison criterion for the solutions 
of two noiseless systems with comparable initial conditions. A necessary con- 
dition for the explosion of the solutions of the problem under investigation is 
provided in this section, together with an illustrative counterexample and a par- 
tial converse. In Section |4l we give upper and lower bounds for the value of the 
time of explosion of our integral equation. Finally, Section [5] provides estimates 
of probabilities associated to the time of explosion of a system in which the 
noise is the absolute value of the Brownian motion. 



2 Osgood's test 

Let M denote the set of extended-real numbers. Throughout, a, b : [0, oo) — > R 
will represent positive, continuous functions, while the function g : [0,oo) — ^ R 
will be continuous and nonnegative. For physical reasons, the function g is called 
a noise. In this work, xq will denote a positive, real number, and X : [0, oo) R 
will be a nonnegative function whose dependency on i > is represented by Xt- 
We are interested in establishing conditions under which the solutions of the 
integral equation 

Xt^xo+ f a{s)b{Xs)ds + g{t), t>0, (2) 
Jo 
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explode in finite time. More precisely, we define the time of explosion of X as 
the nonnegative, extended-real number = sup{t > : Xt < oo}. In this 
manuscript, we investigate conditions under which the time of explosion of X 
is a real number. 

Letting Yt = Xt — git), one sees immediately that the problem under con- 
sideration is equivalent to finding the time of explosion of the solution Y of the 
equation 

Yt = XQ+ [ a{s)b{Ys + g{s))ds, t>0. (3) 

As a matter of fact, = Tj . From this point on, this common, extended-real 
number will be denoted simply by Tg for the sake of briefness. 

Remark 1 It is worth noticing that ([3]) can be presented in differential form 
as the equivalent, initial- value problem 



dt 

Yo = xo 



-^^ait)biYt+g{t)), i>0, 



a problem for which the existence of solutions is guaranteed, for instance, when 
b is locally Lipschitzian and a is regulated (see (10.4.6) in 

Let r be a real number such that < r < xq. We define the functions 
^ : [0, oo) ^ M and B-r : [xq - r, oo) M by 



Ait) 



I ais)ds and BJx) = / tt^- 

Jo Jxo~r KV 



For the sake of convenience, we let B be the function Bq. Evidently, both of 
these functions are nonnegative, increasing and continuous, and so are their 
inverses. On the other hand, if r > —xq, we let Br : [a;o + r,oo) — M be given 
by Br = B-r- For every fixed x > xq, we define B^ : [—xo,x — xq] — >■ R by 
B^ir) = Brix); we prefer this second notation in either case. Additionally, we 
define /3 : [-a;o,oo) -J> R by Pit) = B°°it). All of these functions and their 
inverses are nonnegative, continuous and decreasing in their domains. 

Lemma 2 (Generalized Osgood's test) The initial-value problem 

I ^=^myit)), t>o, 

I 2/(0) = Xo, 

has a unique solution given by y(t) = B^^iAit)), for t < j4^^(_B(oo)). The 
solution explodes in finite time if and only i/ i?(oo) < j4(oo), in which case, 

Ty = A-\Bi^)). 
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Proof. The function y{t) = B-'^{A{t)) is evidently a solution of Ad- 
ditionally, expressing the differential equation in ([5]) as y' {s)/b{y{s)) = a{s), 
integrating both sides over [0,t] and performing a suitable substitution, we ob- 
tain that B{y{t)) = A(t), whence the uniqueness follows. Moreover, y{t) is real 
if and only \it < A'~^{B{oo)). 

Now, if the solution of ([S]) explodes at the time t,, < oo, then B{oo) = 
A{te) < A{oo). Conversely, the number A^^{B{oo)) is real, so that 

B{y{A'\Biocm - AiA-\Bioo))) = B{^). 

This implies that T-f < A^^{B (oo)), and the opposite inequality follows from 
the fact that the solution of (O exists for t < A~^{B{ooj). ■ 
As a consequence, the solution of ([5]) is nonnegative, continuous and in- 
creasing on [0,TJ'), and so is its inverse on [a;o,oo). Likewise, the function 

B : [xo, oo) — > R, given by the formula 

— ds 

^^^^ " 4 bis + giY-^isW 

is nonnegative, continuous and increasing. 

Corollary 3 The solution of ^ can be expressed asYt = B ^^[A(t)), for every 
t < A~'^(B{oo)). 

Proof. The proof runs as in Osgood's test. ■ 



3 A comparison theorem 

Theorem 4 (Comparison criterion) Let < xq < xi, let b be non-decrea- 
sing, and assume that the functions u, v : [0, oo) M satisfy 

v{t) > xi + a{s)b{v{s))ds and u{t) — xq + I a{s)b{u{s))ds, t > 0. 
Jo Jo 

Then, v{t) > u{t) for every t>0, and T^ < A-'^{B (oo)). 

Proof. It is sufficient to show that v > u because, in such case, T^ < T" = 
A~^{B{oo)). Assume first that xq < xi, and let iV = > : u{s) < v{s),s G 
[0,t]}. The set N is nonempty, so T = supA^ exists in M. If T were a real 
number, then 

L = lim (w(r + e)-'u(T-Ke)) > xi -xo-f lim / a{s) \b{v(s)) - b(u(s))] ds 

e-i.O+ £^0+ Jf 

by the fact that v{s) — u{s) > 0, for every s E [0, T]. It follows that L > xi — xq. 
By definition, there exists S > such that v{T + s) — u{T + s) > for every 
s G [0,(5), whence it follows that T + | G A, a contradiction. Consequently, 
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u{t) < v{t) for every t > 0. Now, in case that xq — xi, the solution of the 
equation 

Ur{t) = Xq — r + / a{s)b{ur{s))ds, < r < xo, 
Jq 

satisfies v{t) > Ur{t), for every t > 0. Using Osgood's test and the continuity of 
B^^, we obtain 

v{t) > lim Ur{t) = lim B^-^iAit)) = B~\A{t)) = u{t). 

r->0+ r-i.O+ 



Theorem 5 Suppose that b is non- decreasing, and B{oo) < A{(X)). Then the 
solution of (^) explodes in finite time. The time of explosion of Y is t,, = 
A-i(B(oo)). 

Proof. The fact that b is non-decreasing yields 

Yt = xo+ f a{s)b(Ys + g{s))ds > xq + f a{s)b(Ys)ds. 
Jo Jo 

Theorem H] gives that Tg < A^^{B{oo)) when we compare Y with the solution 
of Yt ^ xo + J*a{s)b{Ys)ds. On the other hand, A{te) = B{oo) < B{oo) < 
A(oo), which implies that te is real. The expression of te and Corollary [3] yield 
= B{YtJ = A{Te). We conclude that T^ = t^. ■ 
Intuitively, it is not generally true that the explosion of the solutions of ([4]) 
in finite time is a sufficient condition for the inequality B{oo) < A{oo) to be 
satisfied. This assertion follows after noticing that A and B do not depend of 
the noise function; however, the time of explosion does. We will establish our 
claim precisely through the following counter-example. 

Example 6 Let xq = 1, and let a, b and g be given by the expressions a{t) = 
e~*, b{t) ~ jt^, and g{t) = e*, for every t > 0. Expanding the expression 

{Y, + e^f in (0]), we obtain Yf > 1 + j /J Y^ds. Then Ft > (1 - \t)-^, which 
implies that Y explodes in finite time. However, B{oo) = 2 > 1 — A{oo). 

The following result is a partial converse of Theorem \E\ We let g{t) = 
sup{5(s) : s S [0,t]}, for every t>0. 

Proposition 7 Suppose that b is non- decreasing, and that 

/oo 
a{s)ds. 

If the solution Y of ^ explodes in finite time, then B(oo) < A(oo). 
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Proof. By Corollary [31 B{oo) ~ B{Yt^) — A{Te). Since b is non-decreasing 
and g{Y~^) < g{Te) for every s G [xq, oo), we obtain 

< A{T,). 



Separating the integral in the definition of B{oo) as the sum of the integrals over 
the intervals [xq, xq +'g{Te)] and [xq +g{Te), oo) , using the facts that b is positive 
and non-decreasing, and employing the last inequality and the hypothesis, we 
obtain 



-B(oo) < 



XO Kxo) Jxo Ks+g{Te)) 



4 Approximation of the explosion time 

It is important to notice that the time of explosion of F, as given by the The- 
orem [5l presents the disadvantage of depending on the solution Y itself. In 
this section, we will derive some approximations to Tg which do not present 
this shortcoming. For the remainder of this manuscript and for the sake of 
convenience, we let T — A'~^{B{oo)). Throughout this section, b will be a 
non-decreasing function. 

The Comparison criterion and Osgood's test yield that the time of explosion 
of the solution F of ([3]) satisfies Tg < T. On the other hand, 

Yt<xo+ f a{s)b{Ys+g{T))ds, 
Jo 

and the Comparison criterion leads us to conclude that 

(/3(5(T))) < Te < T. (6) 

In general, the function b : [0,oo) — > M is sub -multiplicative if there exists a 
positive constant c such that b{xy) < cb{x)b{y), for every x,y > 0. Evidently, 
exponential and power functions are sub-multiplicative. 

Suppose that 5 is a sub-multiplicative function, and let c be the positive 
number provided by the definition of sub-multiplicativity. In the following, it 
will be convenient to define the function A : [0, 00) — > M by 

A{t) = c [ ais)b ( —gis) + l] ds. 
Jo \xo J 

This function is nonnegative, continuous and increasing and, thus, it is invert- 
ible, and has a continuous and increasing inverse. 
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Proposition 8 Let b be a sub-multiplicative function. Then > A ^{B{oo)). 
Proof. Since Yg > xq for every s > 0, we obtain 



Ys g{s) 

< Y,g{s) (^ + ^\^ Yg (^g{s) + l) . 

\xo g[s)J \xo ) 

Monotonicity and sub-multiplicativity of 5, along with (jS]), yield 
Yt<xo + cJ^ a{s)b + 1^ b{Ys)ds. 

The conclusion of the proposition follows now from the Comparison criterion 
and Osgood's test. ■ 



5 An application 

Throughout this section, we consider the stochastic differential equation ([2]) 
with noise function \Wt\, where W is the Brownian motion. The noise function 
is taken in absolute value in view of physical considerations on the dynamics of 
cracks growth under fatigue loading. In fact, it has been established experimen- 
tally that cracks in the subcritical stage grow with a velocity that increases with 
the crack length [13]. The governing equation is called Paris' equation, and it 
is a power law (which is a sub-multiplicative function) in which the exponent is 
determined empirically. As a matter of fact, it has been established that Paris' 
law is valid for a wide range of materials [TJ [2j [TOl [14] . 

For the remainder of this work, we let represent the probability that a 
random variable with standard normal distribution assumes values in [0,a;], for 
every x > 0. 

Proposition 9 LetO <t <T. Then 

Proof. We use here the first inequality of ([6]). Notice that 

PiTe<t)<P (^A-\P{\Wt\)) <t)=P (\Wt\ > r\Ait))) . 

Equation (8.4) in [8] completes the proof. ■ 
For every nonnegative, real number r, we let Tr = inf{i > : \Wt\ = r}. 
Evidently, < r, for every s G [0,Tr]. 
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Proposition 10 Let <t <T. For every r > 0, 



1 - $ (r/^A-^B{Br\Am)) 

P (Te < t\Tr <T)< ^ ^ ^. (8) 

1 - $ [r/VT] 



Proof. Notice that \Wt\ > r whenever Tr < T. Moreover, Osgood's test and 
the Comparison criterion imply that Yt > B^^{A{t)), for every t > 0. Using 
©, we obtain A{Te) > B^^^^{oo) > Br{oo) > Br{B-^{A{Tr))). Therefore, 

P (Br{B-\A{Tr))) < A{t) 

P (Te < t\Tr <T) < — ^ 



P(Tr < T) 

p{Tr<A-\B{B-\A{t)))) 

P{Tr < T) ■ 

The conclusion follows now from [5] as in Proposition [TUl ■ 
On physical grounds, the function Y may represent the temporal behavior of 
the transversal length of a crack failure on some material [T^. In this context, 
the parameter represents the initial, transversal length of the crack, and L is 
the transversal length of the material. For practical purposes, one may think of 
the wing of an airplane which has a fixed transversal length, on which there is 
a crack with known initial length. In such case, one investigates the dynamics 
of the length of the crack with respect to time, in order to conduct preventive 
maintenance on the wing and avoid possible accidents [3]. 

Proposition 11 If L > xq, then 

P {Y-' < t) < I ~ <^ ( ^^''^^^^^ 



T J 

Proof. Let Y the solution of Yt — xq + Jq a{s)b{Ys + \ WT\)ds, for every < 
t <T. By Osgood's test and the Comparison criterion, B^ (^(0) —Yt> Yt. 

\Wt\ 

Once again, the conclusion is reached using [5] in the right-most end of the chain 
of identities and inequalities 

P{Y^'<t) < p{Y^'<t) 

= P {b^w^\{L) < A{t) 
= l~p(\WT\<B-^lt^{L) 
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Example 12 Let xq, uq and a be positive numbers, and let a{t) = qq and 
b{t) — for every i > 0. Observe that A{t) = aot and, for every r E [0,xo] 

and every x > xq, 



Br{x) = - 



SO that T = {aaox^y^. By 0, 



1 1 

(.To — r)" a;" 

(aaoi)-!/" - 



xa 



Vt 



(9) 



for every < t < T. In order to estimate the value of t for which < t with a 
probability of at most 0.05, Equation © yields 



1 - $ 



(aaot) 



-l/a 



Xq 



Vt 



< 0.05 



whence it follows that t = -^[xq + \/T^ ^(0.95)] Proposition |8] and mono- 
tonicity on the integrand imply that 



P{Te<t) < P (^B{oo) < A{t) 

< p(^<i^ao{^m\ + i 



l+a 



ds 



< p 



1 



1 _ $ 



<aot[ — 11^*1 + 1 



-1 



This last estimate of P(Tp < t) is better than that given by (jH]), in view of the 
fact that aaQX^t > 1. 
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